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Abstract. A CCD photometric study of the dense galactic open cluster M37 is presented and discussed. The 
majority of the analysed data are time-series measurements obtained through an Re filter. The observations 
were carried out on seven nights between December 1999 and February 2000, and have led to the discovery of 
7 new variable stars in the field. Three of them have been unambiguously identified as W UMa-type eclipsing 
binaries, while two more are monoperiodic pulsating stars, most probably high-amplitude 8 Scuti-type variables. 
The remaining two stars seem to be long-period eclipsing binaries without firm period determination. Johnson B 
and V frames have been used to construct a new colour-magnitude (CM) diagram of the cluster, and to find the 
locations of the new variable stars. The pulsating variables are most likely background objects. The CM diagram 
is fitted with recent isochrones yielding the main parameters of the cluster. 
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1. Introduction 

Variable stars in open clusters provide an important tool 
to test theoretical predictions concerning stellar struc- 
ture and evolution, because the same distance, age, ini- 
tial chemical abundance and interstellar reddening can be 
assumed for all the stars in the same cluster. Moreover, 
studying cluster variables is one of the most effective and 
productive uses of small and moderate size telescopes. 
Simultaneous photometry of hundreds of stars enables 
very accurate light curves to be obtained. In the recent 
years, substantial development has been made in this area, 
owing to the importance of, for instance, studying stel- 
lar pulsation (Frandsen & Arentoft 1998) or evolution of 
close binaries (Kaluzny & Rucinski 1993) in open clusters. 
Evolutionary effects make different types of variable stars 
appear in specific range of cluster ages, therefore, it is 
highly desirable to make observations of as many clusters 
with different ages as possible. The function of this paper 
is to contribute to this issue with new observations of the 
intermediate-age open cluster M37. 

M37 (=NGC 2099=C0549+325, a 20 oo 
05 h 52 m 3, <5 2 ooo = +32°33') is a very rich open clus- 
ter located in the low galactic latitude (b=+3?09) 
region in the constellation Auriga. Although it is among 



the brightest open clusters in the northern sky, the 
observational material in the literature is quite scanty. 
The last extensive photometry is that of West (1967) 
who presented photographic UBV measurements for 
930 stars in the cluster. The main sequence is rather 
broad and the available proper motions (Jefferys 1962, 
Upgren 1966) do not enable a clear distinction be- 
tween the member and non-member stars. Robin (1982) 
published electronographic plate measurements for 
selected stars in the cluster in the UBV bands for 
calibration purposes. The fifth catalogue of Lynga (1987) 
lists the following cluster parameters: m— M=ll™63, 
E(B-V)=0? 1 31 and log t=8.30. The most recent update 
of these parameters was given by Mermilliod et al. (1996) 
who studied the red giant content of this cluster and 
found m-M=li m 50, E(B-V)=0 rn 29 and log t=8.65. 
In contrast to the above, the WEBDA catalogue 
(http://obswww.unige.ch/webda) lists m— M=ll™67, 
E(B-V)=(T302, log t=8.540, however, the differences of 
the quoted values do not exceed the expected uncertain- 
ties. To our knowledge, there has been no variable star 
survey in the literature despite the relative proximity of 
the cluster. This neglect might be caused by the strong 
concentration of stars that makes it difficult to separate 
individual stars. 
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Table 1. The journal of observations 



Date 


B 


V 


R 


length 


*X3 




frames 


frames 


frames 


(min) 


£ 


1999 Dec. 27/28 


- 


- 


120 


190 


> 


1999 Dec. 29/30 


3 


3 


30 


100 




1999 Dec. 30/31 


- 


- 


187 


280 


(X 


2000 Jan. 28/29 


- 


- 


145 


340 


c_ 


2000 Jan. 30/31 


- 


- 


43 


90 


> 


2000 Jan. 31/Feb. 01 


- 


- 


58 


150 




2000 Feb. 01/02 


- 


- 


75 


170 




Total: 


3 


3 


658 


1320 





2. Observations and data analysis 

CCD BVR photometric observations were carried out on 
7 nights between December 27, 1999 and February 1, 2000 
at the Piszkcstcto Station of the Konkoly Observatory 
(Hungary). The instrument used was the 60/90/180 cm 
Schmidt telescope, equipped with a Photometries AT- 
200 CCD (1536x1024 pixels, KAF-1600 chip with UV- 
coating). The angular resolution was l'.'l /pixel, yielding 
a field of view of 29' x 18'. Time-series observations were 
made through an R filter (60 and 120 s exposures), while 
three B and V filtered frames were also taken on December 
29, 2000 in order to obtain colour information on the de- 
tected stars (20 s, 60 s and 300 s exposures). The choice of 
using an R filter in the time-series observations was due to 
an unfortunate technical difficulty, that made other filters 
unavailable on other nights. The full observing log is given 
in Table 1. 

The basic image reduction was carried out using the 
standard tasks in IRAF. The flat field corrections made 
use of sky-flat images taken during the evening twilight. 
We performed psf-photomctry with the daophot package 
of IRAF using the Moffat point-spread function, which 
turned out to be well suited for fitting the slightly dis- 
torted stellar profiles (caused by the optical layout of the 
Schmidt system). The magnitudes were calculated rela- 
tive to the ensemble mean of a selected set of stars, with 
an arbitrarily chosen zero-point. Typically, we detected 
about 3000 stars on the R frames. The co-ordinate system 
of each frame was corrected for two-dimensional shift, ro- 
tation and optical distortions, relative to a master frame. 
The master frame is the one which contained the largest 
number of stars, and is likely to contain every stars de- 
tected in the series. Although there is still some small 
chance of loosing a few stars that were undetected on the 
given starting image, the overwhelming majority of the 
stars are properly identified. We could identify 2530 stars 
on at least 20 frames. The last step in data filtering was 
to reject those stars that were closer than 30 pixels to the 
image edges. The finally adopted data set contains 2323 in- 
dividual light curves. We note, that R data on Dec. 29/30, 
1999 (30 points) were excluded from the photometric anal- 
ysis because of their larger noise level. 
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Fig. 1. A comparison of the V magnitudes and B— V 
colours obtained in this study with data by West (1967). 
The diagrams are based on 580 common stars. The solid 
lines denote equality. 



The weather conditions were not photometric when 
the three-colour observations were obtained, therefore, we 
could not perform standard all-sky photometry in order to 
standardise the observed colours and magnitudes. Thus, 
we have tied our data to the photoelectric UBV pho- 
tometric secondary standards in the cluster taken from 
Table 2 in West (1967). The colour-dependent transfor- 
mation coefficients e and /i (AV = Av + eA(B — V), 
A(B — V) = /j,A(b — v)) were determined by observing the 
standard sequence of Schild (1983) in the open cluster M67 
(e=0.10±0.03, n= 0.83±0.03). This \x value indicates fil- 
ters or a detector quite far from the standard, fortunately, 
the standard star observations gave well-defined linear 
transformation equations with no indication of second- 
order terms. We have compared the resulting V and B— V 
data with those of West (1967) and can unambiguously 
identify 580 stars common to both sets which are plot- 
ted in Fig. 1. There are some slight systematic differences 
between the V data which are characterized by the mean 
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Fig. 2. The scatter of individual stars as a function of 
V magnitude. Note the absence of the brightest stars (be- 
tween V=ll™0...12? 1 0) in the cluster field, which were sat- 
urated in the time series data observations. 



deviation of V(p.p.)-V(W67) and the rms around that 
mean in certain magnitude ranges. The results are: for 
V= 10-12 mag the mean deviation is — 0™032 with an rms 
(F043; for V= 12-14 mag +CP044, m 049; for V=14 16 
mag +0 m 015, m 067; for V= 16-18 mag -0 m 013, 0™098. 
The colour data are also affected by these systematic dif- 
ferences. The largest deviations occur for the reddest stars, 
our values arc slightly bluer than those of West (1967). 
This inconsistency is due to the fact that red stars be- 
yond B— V«l™50 were essentially extrapolated with the 
standard transformations. The two reddest photoelectric 
standards in Table 2 of West (1967) have B-V l m 626 and 
1™559 and such red stars are very likely variable stars 
(Jorissen et al. 1997). In summary, the uncertainty of 
standard V and B— V data changes from 0™03 to 0™10, 
depending on the colour and apparent brightness of the 
stars. Since we could not standardise the R-band data, 
their actual values contain the arbitrary zero-point. 

The search for variable stars was performed by check- 
ing the light curve statistics and by visual inspection of 
the suspected variable time-series data. We have calcu- 
lated the scatter of individual light curves for every star 
(in sense of (m b s — (mobs))) an d it is plotted as a function 
of V magnitude in Fig. 2. The identified new variables are 
also labelled. The main ridge helped us to characterise the 
accuracy of the brightness measurements, which changed 
between 0™01-0™06 over a range of ^6 magnitudes. As 
can be seen in Fig. 2, besides the strongly deviant points, 
there is a group of slightly deviant points over the main 
concentration. We have checked all of them and found that 
those stars are in the most crowded fields in the cluster, 
where the applied psf-photometry may suffer from larger 
photometric errors. That is why all of the data with sug- 
gested larger scatter was inspected visually. Seven stars 
remained as unambiguously variable stars and these are 
identified in the 2x2' finding charts in Fig. 3. 



The light curves of new variable stars were further 
analysed to search for possible periodicities. For this, 
we performed conventional Fourier analysis of the data 
with Period98 (Sperl 1998). However, in case of eclips- 
ing binaries, the phase dispersion minimisation (PDM - 
Stcllingwcrf 1978) is more suitable for inferring periodic- 
ities. Thus, PDM spectra enabled an independent check 
of the results. Original data are available electronically at 
CDS-Strasbourg. 



3. Results 

The analysis consisted of the following steps: i) construc- 
tion of the colour-magnitude (CM) diagram and physi- 
cal parameter estimation of the cluster; ii) variable star 
classification based on the light curve shapes; in) period 
determination for stars with well-sampled light curves; iv) 
membership investigation (considering the location of new 
variables on the CM diagram and using various period- 
colour-luminosity relations that are valid either for W 
UMa or 5 Scuti stars). 

3.1. The colour-magnitude diagram 

The physical parameters of the cluster have been esti- 
mated by a isochrone fitting of the CM diagram. Since 
we could not determine the interstellar reddening inde- 
pendently this parameter has been also included in the 
fitting procedure. The isochrones were taken from Bcrtelli 
et al. (1994) and the initial parameters were adopted from 
Mermilliod et al. (1996). The result is shown in Fig. 4, 
where the CM diagram is plotted with the two closest 
isochrones z=0.02). Note the strong field star contami- 
nation, which is caused by the position of the cluster in 
the galactic plane. The overall morphology is well repro- 
duced, except that the models are too red for the red gi- 
ant stars. Note that this diagram is essentially the same 
as Fig. 11 in Mermilliod et al. (1996), the only difference 
is that our graph is based on a much larger data set which 
extends to ^2 mag below the faint limit in Mermilliod et 
al. (1996). Therefore, the resulting parameters are very 
similar to those of Mermilliod et al. (1996): the colour 
excess is E(B-V)=(F29±(F03, the distance modulus is 
m-M=ll m 48±0 m 13 and log £=8.65. These values were 
adopted in this paper. 

3.2. The new variable stars 

In the following section we describe the new variables 
and their light variations. Eclipsing binaries and pulsating 
variables are discussed separately. Thanks to the moder- 
ately good phase coverage, five of the seven new variables 
have unambiguous periods and classification. The remain- 
ing two stars have characteristic light variation of long- 
period detached eclipsing binaries, i.e. long constancy in- 
terrupted by short fadings. We have given the stars the 
VI... V7 designations, where the numbers increase with 
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Fig. 3. Finding charts for the newly discovered variable stars. Each frames are 2x2' wide, north is up, east is to the 
left. 



the average magnitude. Their basic data are summarised 
in Table 2. 



0.5585±0.0001 days. The phase diagram is presented in 
Fig. 6. 



3.2.1. Eclipsing binaries 

VI: There was only one fading of m 24 on JD 2451572. 
Unfortunately, we did not observe the subsequent bright- 
ening, thus no epoch of minimum could be determined. 
The first night of observations revealed an ambiguous 
brightening, which may be the ascending branch of a pri- 
mary minimum. The light curves are plotted in panel (a) 
of Fig. 5. The larger scatter in subpanel a3 was caused by 
the bad seeing which increased the light contamination 
from the close bright neighbours (see its finding chart in 
Fig. 3). 

V2: We observed a 0™23 deep minimum on JD 2451540 
which occured at HJD(min)=2451540.518. The other ob- 
servations showed no significant variation, thus no period 
could be determined. The light curves are shown in panel 
(b) of Fig. 5. 

V3: The amplitude is 0™32, while the inferred period is 
0.4224±0.0001 days. The light curve shape is very char- 
acteristic for a contact binary. The phase diagram of the 
whole light curve is shown in Fig. 6. We note the slightly 
different maxima that might be a hint of a spotted stellar 
surface. 

V4: The light curve is typical of an EW-type star with 
an amplitude of O^S. There are interesting humps on 
the light curve suggesting again the possibility of pho- 
tometric surface inhomogeneities. The inferred period is 



V7: This is the faintest variable star detected by us in the 
cluster, and that is why its light curve contains only 185 
points (it has been lost from many frames by the star de- 
tection algorithm). The observed amplitude is 0™55, while 
the period is about 0.35 79 ±0.0001 days. The phase dia- 
gram shows a typical W UMa-type variation (Fig. 6). 

3.2.2. Pulsating variables 

V5: The short-period asymmetric light curve with an am- 
plitude of 0™32 is typical of pulsating stars. The period is 
0.2800±0.0005 days. Since it is located very close to the 
edge of the field, only 5 nights of data were analysed and 
their phase diagram is shown in Fig. 6. 



V6: The light variation is that of a monopcriodic 
high-amplitude S Scuti-star (P=0.109835±0.000005 days, 
A=0?M5). The phase diagram (Fig. 6) also suggests a re- 
markably repetitive light change. 

We have also studied the possible multiperiodic na- 
ture of V5 and V6 by a standard Fourier analysis of their 
light curves. After prewhitcning with the main frequen- 
cies (and their harmonics), the residual frequency spectra 
did not show any significant peak. However, because of 
the relatively short span of observations we cannot firmly 
conclude the strict monopcriodicity of these variable stars. 
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Fig. 4. The colour-magnitude diagram of M37 with two isochrones (solar composition) taken from Bertelli et al. (1994). 
The morphology is well reproduced, although the red giant models are slightly too red. The positions of new variable 
stars are marked with the up triangles. The brightest one is VI (type EA:), the next one is V2 (EA:), while the four 
variables at Vw 16 mag are V4 (EW), V3 (EW), V5 (DSCT/RRc) and V6 (DSCT) starting from the bluest one. V7 
is not presented in this diagram because it was not detected on the B frames. 



3.3. Cluster membership of the variable stars 

We have marked the locations of the variable stars on 
the CM diagram in Fig. 4 and must stress that these po- 
sitions, based on a single-epoch observation, are uncer- 
tain to ~ ±0™2, as both the B— V colour and V mag- 
nitude change over the time. This effect is larger for 



the pulsating variables, which usually have more vari- 
able colour indices than eclipsing stars. The locations of 
the eclipsing binaries are consistent with the stars be- 
ing members of the cluster, as they scatter around the 
main sequence. This can be further studied for V3 and 
V4 with the pcriod-colour-luminosity relation published 
by Rucinski & Duerbeck (1997). Their Eq. (5a)(M v = 
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Fig. 5. The light curves for two new suspected Algol-type eclipsing binaries. 
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Fig. 6. The phase diagrams for V3..V7. 



-4.441og P + 3.02(B-y)o + 0.12) predicts M v (V3)=3 m 62 timated, e.g., with help of the period-colour diagram of W 

and My(V4)=3™00 which is about 1-1.5 mag brighter UMa-type systems (Fig. 2 in Rucinski & Ducrbcck 1997). 

than the calculated absolute magnitudes using the clus- However, the above quoted difference does not exlude the 

ter's distance modulus. Note, that when calculating the membership of V3 and V4, because even the calibrating 
absolute magnitude, the unreddened colour should be es- 
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Table 2. The basic data of new variables 



Star 


02000 


(52000 


(V) 


{B- 


V) 


P(d) 


A (mag) 


Epoch(s) 




Type 


VI 


05 h 52 m 20:42 


+32°33'19'.'5 


13 m 78 


0T61 




- 


- 


- 




EA: 


V2 


05 h 52 m 16*60 


+32°28'14'/9 


14 m 98 


0^69 




- 


- 


2451540.518 




EA: 


V3 


05 h 52 m 3£03 


+32°32'41'.'7 


16™07 


(T90 




0.4224(1) 


0.31 


2451575.5083 


(ID 


EW 


V4 


05 h 52 m 5£26 


+32°33'01'.'2 


16 m 17 


(F87 




0.5585(1) 


0.33 


2451576.456: 


(I) 


EW 


V5 


05 h 53 m 00!63 


+32°24'50"8 


16 m ll 


(T96 




0.2800(5) 


0.32 


2451576.500: 




DSCT 


V6 


O5 h 51 m 50!55 


+32°32'34'.'4 


16 m 05 


l m 02 




0.109835(5) 


0.45 


2451540.5367 




DSCT 


V7 


05 h 52 m 39>?32 


+32° 36' 30'.' 9 


17™ 83 


- 




0.3579(1) 


0.55 


2451574.262: 


(I) 


EW 



data set of Rucinski & Duerbeck (1997) contained stars 
with similarly large deviations. 

Contrary to the eclipsing binaries, the pulsating vari- 
ables are too red and too faint to be cluster members. The 
S Scuti instability strip of the H-R diagram (Breger 2000, 
Rodriguez & Breger 2001) is in a bluer and much brighter 
region than that containing V5 and V6 (e.g. the apparent 
brightness range of the instability strip of M37 is about 
ll m 5-14 m 0). 

Petersen & Christensen-Dalsgaard (1999) discussed 
the Hipparcos period-luminosity relation for high- 
amplitude S Scuti stars and concluded that quite good 
absolute magnitudes (±0™1) can be inferred from the re- 
sulting PL-relation. Their Eq. (4) (My = -3.751og P - 
1.969) predicts M v (V5)=0 m l and M v (V6)=l m 6. Both 
values suggest the stars are much further away than M37. 
Unfortunately, farther means redder and the correspond- 
ing higher colour excesses cannot be properly estimated. 

Furthermore, the period and light curve shape of V5 
also suggests the possibility of being a field RRc star in- 
stead of high-amplitude 5 Scuti star. In that case its ab- 
solute magnitude should be about 0™5 (log L « 1.7, e.g. 
Kollath et al. 2000), which yields to a distance modu- 
lus m— M(V5)=15™5. The RRc assumption gives a mean 
(B-V)o colour of 0^2-0^4 (e.g. Lee & Carney 1999) re- 
sulting in E(B— V)«0™7. On the other hand, there are a 
few longer period S Scuti stars known with periods be- 
tween 0.25 and 0.3 days (Rodriguez & Breger 2001), that 
could be evolved massive 5 Scuti variables. Thus, V5 can 
be either a new field RRc star or a high-amplitude, longer 
period S Scuti star. In any case, neither V5 nor V6 seems 
to be physical a member of M37. 

4. Summary 

We have presented the first variable star survey of the rich 
open cluster and have discovered seven new variable stars, 
two pulsating and five eclipsing. Epochs and periods are 
determined for three W UMa-type systems and two pul- 
sating stars. The eclipsing binaries are probably cluster 
members, while the pulsating stars are background ob- 
jects lying much further away than M37. V5 might be 
cither a high-amplitude 5 Scuti or a field RRc star, while 
V6 is classified unambiguously as high-amplitude S Scuti 
star. More observations are needed to determine periods 
for VI and V2, while it is an interesting question whether 



the observed light curve shapes of the contact binaries are 
stable over a longer period. All of the three W UMa sys- 
tems have light curves with different consecutive maxima 
suggesting surface inhomogeneities that can be caused by 
stellar spot activity. 

In addition to the variable star survey, a new 
colour-magnitude diagram has been constructed which 
is based on about 2300 stars. Since the cluster lies 
very close to the galactic plane, the field star contam- 
ination is fairly strong. We have determined indepen- 
dently the cluster's main physical parameters by fitting 
isochrones taken from Bertclli et al. (1994). The best 
fits yield E(B-V)=0 m 29±0™03, m-M=ll™48±0 m 13 and 
log t=8.65. In more convenient units, the cluster's age is 
about 450T0 6 years and the distance to M37 is about 1300 
pc. These parametes are in good agreement with those 
available in the literature. 
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